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Abstract 

The coronavirus responsible for the severe acute respiratory syndrome (SARS-CoV) contains a small 
envelope protein, E, with putative involvement in host cell apoptosis and virus morphogenesis. It has 
been suggested that E protein can form a membrane destabilizing transmembrane (TM) hairpin, or 
homooligomerize to form a regular TM a-helical bundle. We have shown previously that the topology of 
the a-helical putative TM domain of E protein (ETM), flanked by two lysine residues at C and N termini 
to improve solubility, is consistent with a regular TM a-helix, with orientational parameters in lipid 
bilayers that are consistent with a homopentameric model. Herein, we show that this peptide, 
reconstituted in lipid bilayers, shows sodium conductance. Channel activity is inhibited by the anti- 
influenza drug amantadine, which was found to bind our preparation with moderate affinity. Results 
obtained from single or double mutants indicate that the organization of the transmembrane pore is 
consistent with our previously reported pentameric a-helical bundle model. 
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The causative agent of severe acute respiratory syndrome 
(SARS), SARS-CoV, is a member of the family Corona- 
viridae (Rota et al. 2003), which causes common colds in 
humans and is responsible for serious diseases in other 
animal species. The lipid bilayer that surrounds the 
virions typically embeds three proteins; spike (S), mem¬ 
brane (M), and envelope (E) protein. These proteins are 
studied for their important roles in receptor binding and 
virion budding. 

The E protein of the SARS-CoV virus is a small, 76- 
amino acid, integral membrane protein with one putative 
transmembrane a-helical (TMH) hydrophobic domain, 
20-30 amino acids long, flanked by a short N-terminal 
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region of <10 amino acids and a longer C-terminal tail, 
both more hydrophilic. This protein has been found to be 
critical for viral budding in mouse hepatitis virus (MHV) 
(Fischer et al. 1998). Also, in many other coronaviruses, 
coexpression of M and E proteins, which probably interact 
via their cytoplasmic domains in pre-Golgi compartments 
(Lim and Liu 2001), is necessary for morphogenesis (Bos 
et al. 1996; Corse and Machamer 2001). 

Further functional roles suggested for E protein are in 
promoting apoptosis (An et al. 1999; Chen et al. 2001) 
and alteration of membrane permeability when expressed 
in Escherichia coli and mammalian cells (Liao et al. 2004). 
A somewhat related phenomenon is its cation-selective ion 
channel activity (Wilson et al. 2004) observed in vitro, 
which has been reported for a synthetic, full-length, SARS- 
CoV E protein, and also for its N-terminal (40 residues) 
fragment that encompasses the a-helical transmembrane 
domain (Wilson et al. 2004). The channel activity of E 
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protein from two other coronaviruses has been shown to be 
blocked by hexamethylene amiloride (Wilson et al. 2006), 
and this effect has been correlated with inhibition of repli¬ 
cation of these viruses. This points to a possible pharmaco¬ 
logical target during coronavirus infection, although the 
precise function of these putative channels is not well 
known. 

These results suggest that the E protein in corona- 
viruses, specifically that of SARS-CoV, belong to the 
viroporin class of proteins (Pinto et al. 1992; Ewart et al. 
1996; van Kuppeveld et al. 1997; Bodelon et al. 2002; 
Fischer and Sansom 2002; Gonzalez and Carrasco 2003). 
One of the members of this class, M2 from influenza A 
virus, binds the drug amantadine with high affinity (Hay 
et al. 1985). Amantadine has long been used for the 
treatment of influenza A (Fleming 2001). This drug has 
also been used to abolish ion conductance of hepatitis C 
virus p7 protein (HCV), which was also suggested to be a 
viroporin (Griffin et al. 2003). Indeed, amantadine, in 
combination with other drugs, has shown some degree of 
success in HCV clinical trials (Torre et al. 2001; Zilly 
et al. 2002). By analogy, amantadine treatment may be 
beneficial against SARS-CoV infection. 

The channel properties suggested for SARS-CoV E 
protein may be caused by oligomerization of its trans- 
membrane domain (ETM) to form an a-helical bundle. 
Indeed, we have shown recently that the synthetic trans- 
membrane domain of SARS-CoV E, flanked by two 
lysine residues (K 2 -ETM-K 2 ) to improve solubility, insert 
in lipid bilayers as a regular a-helix (Torres et al. 2006). 
Multiple infrared dichroic data were consistent with one 
of the homopentameric models we obtained computation¬ 
ally using evolutionary conservation data (Torres et al. 
2005). We proposed that this a-helical bundle was 
responsible for the channel activity observed. This pre¬ 
diction, however, was in contrast with a short a-helical 
hairpin topology, suggested using a slightly longer ETM, 
which did not contain flanking lysine residues (Arbely 
et al. 2004; Khattari et al. 2006). 

Herein, we have performed experiments to determine 
if peptide K 2 -ETM-K 2 and several of its mutants show 
channel properties consistent with those observed for the 
full-length and 40-residue N-terminal peptide of SARS- 
CoV E (Wilson et al. 2004). Also, we wanted to test if 
amantadine can block this channel activity, and measure 
its binding affinity. Our results suggest that our modified 
peptide, K 2 -ETM-K 2 , explains the channel properties of 
ETM, and can be used as a model to study the channel 
formed by SARS-CoV E. 

Results and Discussion 

The peptides synthesized for K 2 -ETM-K 2 and its mutants 
are shown in Figure 1. The positions of these mutated 


10 15 20 25 30 35 
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WT KKTGTLIVNSVLLFLAFVVFLLVTLAILTKK 

N15A KKTGTLIVASVLLFLAFVVFLLVTLAILTKK 

N15A,V25F KKTGTLIVASVLLFLAFVFFLLVTLAILTKK 
V25F KKTGTLIVNSVLLFLAFVFFLLVTLAILTKK 

V2 5 F,A3 2 F KKTGTLIVNSVLLFLAAVFFLLVTLFILTKK 
F23A KKTGTLIVNSVLLFLAAVVFLLVTLAILTKK 

Figure 1. Synthetic peptides corresponding to the transmembrane 
a-helical domain of SARS-CoV E (ETM) protein (residues 9-35) and sev¬ 
eral of its mutants, flanked by two lysine residues at the C and N termini. 
The mutations introduced are indicated in bold. The numbering {top) 
corresponds to the SARS-CoV E protein sequence. 


residues in the context of the a-helical pentameric bundle 
model we reported previously (Torres et al. 2006), are 
shown in Figure 2. This latter figure shows that residues 
N15, V25, and A32 face the lumen of the pore, whereas 
F23 is oriented toward lipidic environment. 

SDS electrophoreses for all the peptides shown in 
Figure 1 (not shown) showed a similar pattern of 
mobility, consistent with dimers, trimers, and pentamers. 
No significant amount of monomers (~3.3 kDa) were 
detected, even in the case of double mutations. In the 
latter case, it is likely that native monomer-monomer 
interactions are affected; therefore, the bands found are 
likely to represent nonnative interactions induced by 
SDS. These results were identical to those we published 
for the wild-type (labeled WT in Fig. 1) (Torres et al. 
2005, 2006). 

Conductance experiments 

Figure 3 shows the conductance results for our peptide 
K 2 -ETM-K 2 and its mutants. The I/V plots corresponding 
to these experiments, which were used to estimate the 
conductance, are shown in Figure 4. The wild-type (WT 
in Fig. 1) conductance (—45 pS) is consistent with that 
observed in previous reports for the full-length and N- 
terminal fragment (Wilson et al. 2004). This demonstrates 
that the stretch of amino acids used here is responsible for 
channel activity. Also, we point out that this conductance 
is observed in the presence of the two lysine residues at 
N and C termini, which were introduced to increase sol¬ 
ubility and facilitate purification of the synthetic peptide. 

The single-mutants, N15A and V25F, did not show any 
conductance at any of the voltages tested (see Fig. 3), 
consistent with the orientation of these residues toward 
the center of the pore (see Fig. 2) (Torres et al. 2006). The 
single-mutation F23A showed conductance (114 pS), 
consistent with an orientation of this residue toward the 
lipidic environment (see Fig. 2). 

The double-mutants, V25F-A32F and N15A-V25F, 
showed some conductance, 77 pS and 10 pS, respectively, 
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Figure 2. Slices through the ETM pentameric a-helical bundle (Torres et al. 2005, 2006). The residues mutated in this work are 
indicated inside the circles. 


despite the fact that they contain the mutations V25F and 
N15A. These latter mutations on their own, as shown 
above, completely abolish conductance. The high con¬ 
ductance observed in this mutant could be due to non¬ 
specific destabilization of the membrane. 

To test this hypothesis, we compared the ion selectivity 
of WT and double-mutant V25F-A32F (Fig. 5). This 
figure (panel C) shows an asymmetrical current-voltage 
relationship for the wild-type sequence (WT), which 
indicates that it is more permeable to sodium ions than 
to calcium ions. In contrast, the double mutant shows a 
symmetrical current-voltage relationship which indicates 
no cation specificity. 

The effect of the addition of amantadine on conduc¬ 
tance is shown in Figure 6. Amantadine blocks conduc¬ 
tance in the native sequence (WT) and in the mutant 
F23A, suggesting a native, or close to native, organization 
of the monomers in these two samples. However, con¬ 
ductance in the double mutants (NI5A-V25F and V25A- 
A32F) was not affected by the addition of amantadine, 
which may reflect a significant structural change. Two 
situations are possible: (1) Amantadine still binds but is 
not effective in inhibiting ion flow, or (2) amantadine 
does not bind at all. The latter, as stated above, would 
explain the high conductance of the double mutants 
relative to the single mutants by a destabilization of the 
membrane. 

Surface plasmon resonance 

To determine the affinity of amantadine for the pore 
formed by ETM and its mutants in lipid bilayers, we used 


surface plasmon resonance (SPR). Figure 7 depicts the 
amantadine dose-response curves obtained under equili¬ 
brium conditions for the sequences shown in Figure 1. 
The data were fitted using a model that assumes a one- 
to-one binding of amantadine and channel. 

Figure 7 shows that amantadine binds with relatively 
high affinity to native ETM (—10 ! M). This affinity is 
comparable to that observed for influenza A M2, ~10 -3 - 
1CP 4 M, depending on the strain (Astrahan et al. 2004). 
As predicted from the model in Figure 2, significant 
binding was only observed for F23A ( k D ~14 mM), 
whereas only what appeared to be nonspecific aggrega¬ 
tion was observed for N15A, V25F, and A32F, suggesting 
an alteration of the native oligomeric structure in these 
samples. This is also consistent with our observation that 
conductance cannot be inhibited by amantadine in V25F- 
A23F or N15AV25F (Fig. 6), despite the fact that these 
samples show high conductance. Therefore, we conclude 
that the latter is due to membrane permeabilization in 
this case. 

The observation that mutation at N15 disrupts ion 
channel conductance and amantadine binding is expected, 
attending the conservation of a polar residue (S, Q, or N) 
at this position (Torres et al. 2005) in all coronavirus 
envelope proteins. However, the results obtained with the 
other mutations could not have been anticipated in the 
absence of our model. In summary, our results not only 
confirm that the ion conductance properties of SARS- 
CoV E reside in the transmembrane domain, we also 
show that a peptide flanked by two lysine residues (K 2 - 
ETM-K 2 ) (Torres et al. 2006) for which a homopenta- 
meric a-helical bundle model was proposed, is a good 
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Figure 3. ETM ion channel activity in planar lipid bilayers for peptides 
shown in Figure 1. Openings are deviations from the baseline (short dashed 
line). The trans chamber was earthed and the cis chamber was held at 
various potentials between 0 {bottom) and 80 mV {top), in 20-mV intervals. 


model to study its channel properties. Further, we report 
that this peptide incorporated in lipid bilayers binds 
amantadine with relatively high affinity, at least compa¬ 
rable to influenza A M2. This result may have therapeutic 
implications in the treatment of SARS infection. 


Materials and Methods 

Peptide synthesis 

Peptides corresponding to the transmembrane domain of SARS- 
CoV E. from residues 9-35, with two lysine residues flanking 
both N and C termini, were obtained as described previously 
(Torres et al. 2006) from W.M. Keck Biotechnology Resource 
Center. The peptides were cleaved from the resin with trifluoro- 
acetic acid (TFA) and lyophilized. The lyophilized peptides 
were dissolved in trifluoroethanol (TFE), TFA. and acetonitrile 
(1:1:4, v/v/v; final peptide concentration ca. 5 mg/mL) and 
injected to a 20-mL Jupiter 5 C4-300 column (Phenomenex) 
equilibrated with H 2 0. Peptide elution was achieved with a 
linear gradient to a final solvent composition of 10% H 2 0, 90% 
acetonitrile, using a Waters 600 HPLC system. All solvents 


contained 0.1% (v/v) TFA. The resulting fractions were pooled 
and lyophilized. Peptide purity was confirmed by electrospray 
mass spectrometry (ESI), and the sample did not have deletions 
or truncations, with a purity of ~95%. 


Peptide reconstitution in lipid bilayers 

Approximately 1 mg of lyophilized peptide and 10 mg of 1,2- 
dimyristoyl-sn-glycero-3-phosphocholine (Avanti polar lipids) 
were dissolved in 500 |xL hexafluoro-2-propanol (HFIP). 
The solution was mixed at 37°C for 1 h and HFIP was com¬ 
pletely evaporated. The protein-lipid mixture was dissolved in 
filtered and degassed phosphate-buffered saline (PBS), pH 
6.6 (10 mM Na 2 HP0 4 • NaH 2 P0 4 , 120 mM NaCl, and 2.7 mM 
KC1) to a final lipid concentration of 3 mM. The solution 
was mixed for 20 min, resulting in a suspension of peptide 
reconstituted in liposomes. Proper ETM reconstitution in lipid 
bilayers was confirmed by Fourier Transform InfraRed (FTIR) 
spectroscopy (Torres et al. 2000), with an amide I band centered at 
1655 cnf (not shown) consistent with a-helix. The liposomes 
were resized to a uniform diameter of 100 nm by extrusion 
through an Avanti Mini Extruder with a 100 nm polycarbonate 
membrane. 


Surface plasmon resonance 

Binding measurements were performed on a BIAcore 3000 
system, with LI chip as substrate, which allows liposome 
binding (Cooper et al. 2000). The LI chip was cleaned prior 
to each experiment with 40 mM Triton X-100 (20 p,L at a flow 
rate of 20 |xL/min). The liposomes were immobilized on one of 
the flow cells of the chip (80 p,L at a flow rate of 2 p,L/min). 
Excess unbound liposomes was removed with 10 mM NaOH in 
PBS (50 |xL at a flow rate of 100 |xL/min) followed by PBS flow 
for 4 h at a rate of 5 p,L/min. Amantadine was injected (60 p,L at 
a flow rate of 20 ptL/min) in the concentration range of 0.5- 
20 mM on the membrane immobilized on the chip. At the end of 
each injection, after equilibrium, amantadine was allowed to 
dissociate for 1 min under PBS. Subsequently, the surface was 
regenerated with PBS (30 ptL/min flow rate for 15 min). 

Each experiment included the immobilization of empty lipo¬ 
somes (i.e., those without peptide) in one of the flow cells. 



Voltage (mV) 

Figure 4. Current-voltage plots (I-V plots) are shown in the lower right 
panel. Single-channel conductances (slopes of the linear fit obtained from 
SigmaPlot, version 9.2) were 44.5 pS for WT (•), 114 pS for F23A (■), 
77 pS for V25F-A32F (▼), 10 pS for N15A-V25F (A). The data 
represents the average from 4 to 10 recordings. 
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Figure 5. Comparison of ion selectivity (see Materials and Methods), with 
single-channel currents, between wild-type sequence (A) and V25F-A32F 
double mutant (B). Voltages are from —60 mV ( bottom ) to 60 mV (top), 
in steps of 20 mV. Short dashed lines indicate the closed state. The 
corresponding current-voltage relationship is shown in (O- A line is 
drawn only to guide the eye. Data for wild-type (•) and double-mutant 
V25F-A32F (O) are indicated. Vertical bars represent two standard 
deviations. 


which was used as a reference. The mean RU (resonance units) 
of the last 10-15 sec of amantadine injection on the empty 
liposomes’ cell were subtracted from the mean RU, at the same 
time period, in the cell containing the reconstituted ETM. Due 
to the fact that the on- and off-rate constants of amantadine 
binding to ETM were rapid, all thermodynamic analyses were 
based on equilibrium, steady-state measurements (see below). 

The binding results were represented as dose-response 
curves, in which the final RUeq is plotted against amantadine 
concentration. The dose-response curves were fitted according 
to Equation 1, where the association constant is given 


by K A = 


Ir 

*on 

k 0 ff 


RU eq = 


K a ■ [Am] ■ RUmax 

1 +K a - [Am] 


( 1 ) 


and [Am] is the concentration of amantadine injected, RU max is 
the RU when all the ETM channels are saturated with amanta¬ 
dine, and RU eq is the RU at equilibrium for a given concen¬ 
tration of amantadine (Astrahan et al. 2004). From K A , K D can 
be obtained (A D = K A '). The dissociation constant (K D ) was 
derived using the Biaevaluation 3.0 program (Biacore). The 
goodness of fit was estimated as x 2 values, where a good fit is 
estimated as \ 2 < 10. These values were well below 10 in all 
cases (ETM, 1.98; V25F, 0.799; V25F-A32F, 1.31; N15A, 1.69; 
N15A-V25F, 1.3; F23A, 0.621). 

The calculated RU MAX was normalized by multiplication by a 
constant factor (Astrahan et al. 2004) and the fitting curve program 
was then run once again to produce calibrated x 2 , to compare the 
different ETMs \ 2 - The model used to analyze the surface plasmon 
resonance (SPR) data assumed a 1:1 binding stoichiometry 
between amantadine and ETM (Astrahan et al. 2004). 


Ion conductance 

Synthetic l-palmitoyl-2-oleoyl-sn-glycero-3-phosphotidylethanol- 
amine (POPE), 1 -palmitoyl-2-oleoyl-sn-glycero-3-phosphotidylserine 
(POPS), l-palmitoyl-2-oleoyl-sn-glycero-3-phosphotidylcholine 
(POPC), and dimyristoyl-phosphotidylcholine (DMPC) were 
purchased from Avanti Polar Lipids. Decane (>99%) and (N-[2- 
hydroxyethyl]piperazine-N , -[2-ethanesulfonicacid]) were obtained 
from Sigma-Aldrich, chloroform from Sigma-Aldrich, and ethanol 
(>99%) and sodium chloride (NaCl) were obtained from Fluka. 
Ion channel formation of ETM was measured by (Wilson et al. 
2004) forming planar lipid bilayers across a 250-p.m aperture in 
the wall of a Delrin cup separating the aqueous solutions in the cis 
and trans chambers. The working volume was of 1 mL in each 
chamber. A lipid mixture of 5:3:2, POPE:POPS:POPC (50 mg/mL) 
was dissolved in chloroform and dried under N 2 gas. The lipid 
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Figure 6. Effect of amantadine addition (100 (jlM, see arrow) on the 
conductance of the peptides shown in Figure 1, except N15A and V25F, 
which did not show conductance (Fig. 3). 
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Figure 7. Representative amantadine SPR dose-response curves. The data is shown in white circles, and were fitted (solid line) using 
Formula 1. The determined standard deviation of the k 0 did not exceed 20% of the average value in each case. 


cake was then resuspended in decane. A small amount of lipid 
mixture was painted across the aperture of the Delrin cup. About 
14-20 |xg of SARS-CoV ETM protein was added to the cis 
chamber under continuous stirring, until ion channel activity was 
detected, typically after 15-30 min. The chambers contained 500 
rnM NaCl, 5 mM HEPES, pH 7.2 (cis), and 5 rnM NaCl and 5 mM 
of HEPES, pH 7.2 ( trans ). Recordings were performed in asym¬ 
metric ionic conditions using a different buffer. The cis and trans 
solutions were connected via an Ag/AgCl electrode and 1 M 
NaCl, 2% agar bridge. Currents were amplified using a Warner 
Instrument amplifier and digitized at 5 kHz and filtered at 10 kHz 
using pCLAMP 9.2 software (Axon Instruments, Inc.). The trans 
chamber was set as reference and the cis chamber was held at 
different potentials ranging from 0 to 80 mV, in 20 mV incre¬ 
ments. Stirring was stopped after detection of channel activity. 
The ion selectivity experiment was performed with 500 mM NaCl, 
5 mM HEPES (pH 7.2) in the cis chamber, and 500 mM CaCl 2 , 
5 mM HEPES (pH 7.2) in the trans chamber. 
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